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Note on the Lunar Theory. By Professor Ernest 
W. Brown, M.A. 

I have determined the following coefficients in longitude de¬ 
pending on the eccentricity of the Sun. The notation used is 
Delaunay’s. The coefficients are the parts of the form e' f (m), 
and are accurate to the last place of decimals given. 

V = ... — 659"-2375 sin Z'+ 152"*0828 sin (2D — V) — 2i /, *5942 sin (2D + Z') 

+ 1" 2550 sin (4D — l') — o"-i8oo sin (4D + l’) 

+ o"oio5 sin (6D— V) — o" ooi5 sin (6D + Z') 

+ o"-oooi sin (8D — V). 

These agree well with Delaunay’s values when estimate is 
made for the omitted portions of the series in powers of m, the 
greatest difference being o"'o6. 

I hope shortly to publish the coefficients of the form ee } f (m), 
which are amongst the least accurate of all Delaunay’s coef¬ 
ficients owing to the slow convergence of the series in powers 
of m. 

Haverford College : 

1894 May 3. 


On some Points connected with the Integration of the Differential 
Equations of the Relative Motions of Material Systems. By 
E. J. Stone, M.A., E.B.S., Badcliffe Observer, Oxford. . 

In problems of the relative motions of material points subject 
to their mutual gravitational attractions, when the masses are 
denoted by 

1 : ^ : V- • • • v v > 

and the corresponding coordinates, referred to the unit-mass as 
origin, by 

£l> ^l> Cl > £2* £> j &c., 


the three differential equations for the motion of each mass assume 
the well-known forms— 


(0 


5f=-/(«+»’.) J+J/-. 

§=-/(i + ^ + j/. 


dR 
<*£1 
dR 
dr) 1 
dR 


&c. = &c. 
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where p x 2 = ii 2 + 1 2 + C 2 ; and /denotes the accelerating effect 
of the unit-mass at the unit-distance in the unit-time. 

The symbol / denotes summation extending over all the 
material points of the system. If any of the masses are so small 
that their accelerating effects are insensible, the corresponding 
symbols v may be neglected in the summation ; but if included 
in the mathematical work, numerical values, sensibly zero, will 
be found for them from a discussion of the observations when 
accuracy has been secured. 

If found convenient the relative motions of any two material 
points about the unit-mass may be resolved, in the mathematical 
treatment, into that of the motion of their common centre of 
gravity about the unit-mass, and their motion about the common 
centre of gravity ; and any other exact algebraical transforma¬ 
tions of the equations may be introduced and used without 
modification of the general principles upon which the solution of 
such physical problems depend. 

The 377 differential equations for the motion of such a system 
contain (jp -f 1) constants / v u v 2 . . . r p ; and their general 
integrals, if determinable, will contain, in addition to these, the 
6 p constants which are required to give the necessary generality 
to the results, so that they may include all the possible cases of 
the relative motions of such a system of material points. If the 
adopted units of length and of time are independent of the 
particular system of relative motions under discussion, the whole 
of the (7 p + 1) constants will have to be found from observations 
of that system, and their determination will be sufficient for the 
complete solution of the physical problem. But if the times are 
found from the relative motions themselves, the determination 
of the “ time ” from observation will replace the determination of 
the two constants, to which definite values are assigned in the 
equations to render the variable t determinate. 

If we express the coordinates under the forms 


and thus introduce six functions to express the position of each 
mass relatively to the unit-mass, we shall require for each mass, 
besides the three differential equations of motion, three additional 
equations, differential or otherwise, to render the coordinates 

x + 5.r, y -f Zz 

determinate. When additional functions, such as ex, ly , cz, 
are introduced, the functions x, y , z may be called primary 
functions, and the functions lx, ly, Iz, auxiliary functions. When 
the primary functions x, y , z are known, the auxiliary functions 
will be determinable from the differential equations. 


(2) • 


!*!!£= —/■(,+»)—— 
dp 1 K p 3 



V 


* d(x + 8x) 


&C. &(\ 
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But as no loss of generality necessarily results from the intro¬ 
duction of these auxiliary functions, the general integrals found 
by adding them to the primary functions will include the rela¬ 
tive motions of any system of (p + i) material points however 
distributed, and moving at a particular instant; and from the 
equal generality of the results, the real difficulties of the general 
integration of the equations (2) with selected primary functions 
will not be less than those which would be experienced in the 
integration of the original differential equations of motion. But 
if we have found in any ivay values of the functions x ly y j, z u 
x. n V 2> #2? &c., which represent approximately the motions in 
some particular case of mass distribution and motion with which 
we are concerned \ and we can, therefore, rightly assume that in 
this particular case the ratios 

Sx 82/ 5 z 

—» 1 > 

r r r 

where 

r = Vd'+y* + 

are always small fractions; we can, in finding the auxiliary func¬ 
tions required for the complete integration of the differential 
equations of motion for the particular system with which we are 
concerned, restrict our investigations to those cases in which the 
ratios 

82c 8y 80 

■ > » 

r r r 

are small fractions ; and the integration of the differential equa¬ 
tions (2) can then be effected by successive approximations in 
which different orders of the fractions 

8x Sy hz 

9 9 

r r r 

are neglected in succession. 

We shall not obtain by following this method the general 
integrals for the relative motions of any system of (;; + i) 
material points, but those for the special system with which we 
are practically concerned. 

If the primary functions x , y , 2 adopted for any moving mass, 
V, are definite functions of t , and contain no arbitrary constants 
to be found from observation, the six arbitrary constants 
required to give the necessary generality to the results must 
be introduced in the integration of the differential equations (2) ; 
and in this case the general integrals of the auxiliary functions 
Sx, Sy , Ss will have to be found. 

But if the necessary constants have been introduced in the 
primary functions x, y , z, it will only be necessary to find par¬ 
ticular integrals of the differential equations (2), and the work 
will be greatly facilitated in those cases in which the ratios 

Sx 8y 8 z 

9 - 9 * 

r r r 
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are small fractions ; but in approximating to the expressions for 
these particular integrals it will be impossible to assign to the 
arbitrary constants contained in the functions x, y, z approximate 
numerical values, unless the errors which arise from the use of 
these approximate values are at every stage of the mathematical 
work cf the order of the neglected terms 


dx 

— 5 

r 


By Bz 

“ ^ J '— • 

r r 


if we adopt, as is usual in the problems of gravitational astro¬ 
nomy, the complete integrals of the differential equations 


( 3 ) 


d 2 x 

dt l 


= -/(!+*’) 


d pL=-f{ I4f) 


d?z 

~dd‘ 


-/(!+") 



for the primary functions, we can express x, y , % under the. 
forms 

x = r[c as v . cos Q — sin v . sin £3 . cos 7] 
y =r[cos v .sin £3 + sin v . cos S3 . cos 7] 
z=r . sin v . sin 7 

r = a(i — c cos <p) : tan Mil = —~- t , tan (p 

<p-e sin <p= III) . (2 + t) 

where a , e, o, r, p, y are six constants introduced in the integra¬ 
tions. 

But the primary functions x, y, z thus obtained are not the 
integrals of the differential equations of motion unless the ac¬ 
celerating effects of the masses v l5 i/ 2 , v- 6 , . . . v v are insensible, 
and this is not the case in the problems with which we have to 
deal in astronomy. It will not, therefore, be possible to obtain 
•with sensible accuracy the numerical values of the constants- 
contained in the differential equations of motion and their in¬ 
tegrals by equating the functions x, y, z with the coordinates of 
the moving masses as found from observation. If this is 
attempted, and definite values are assigned to x, y , % in the 
formation of the differential equations (2), the general integrals, 
instead of the partial integrals, of these equations will be required 
to ensure accuracy. If, however, we accept the primary elliptic- 
integrals for the functions x, y, z r and use them, in their general 
form , in the formation of the differential equations 

+ Br .r 1 r , <rZR 

P :{ r J J + J ' ' V * d{x T£) 

&c. = &c. 


( 4 ) 


(l-ox 

llF 


= -/(I+*0 
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the partial integrals Sx, Sy, Sz found from these equations, free 
from any additional arbitrary constants, will by their addition to 
the primary functions x , y, z give solutions of the differential 
equations of motion. But the integration of the equations (4) 
will not lead to the required integrals of the differential equations 
of motion unless the functions x, y , s which are actually used for 
the formation of the differential equations (2) are integrals of 
the primary differential equations (3). This point is of more 
practical importance than may at first sight appear to be the 
case, and it is liable to be overlooked in carrying out the 
laborious operations required for the determination of the 
auxiliary functions Sx, Sy, Sz for the different masses. The same 
units of length and of time will have to be adopted in the forma¬ 
tion of the differential equations of motion for all the material 
points whose relative motions are under discussion ; or, what 
practically comes to the same thing, if different units are employed 
the ratios of the different units must be accurately known, and 
allowed for in the formation of the differential equations. But 
when the same units of length and of time are adopted, the- 
numerical value of the factor f which appears in each of the 
primary functions, will require that the same numerical value 
should be assigned to it in all the primary functions and re¬ 
tained throughout the investigations of the auxiliary functions. 
Sx, Sy, Sz. When, therefore, angular velocities 

. . . n p 

are introduced to replace the complex functions 

//(I +vf) / f{l+v p ) 

\f a x s ’ \/ a* ’ \j a p 3 


the primary differential equations will not be satisfied unless we- 
have identically for each material point 

jy _ d-^n-^ . j* _ d^n^ . » _ d^pivp, 

I + 1+^2 1 + v p 

and, in order that the units of length and of time may be the 
same throughout our mathematical investigations of the motion 
of the system, the same numerical value must be assigned to f 
in all the primary integrals, and we must have, therefore, the 
following identities : 

#] 3 W, 2 _ d^lf __ dp S 7lp 2 

I + y, I + V 2 I + Vp 


In the equations (6) a x , a 2 , ... and n u n . 2 , ... denote constants 
introduced in the integrations which have to be found from 
observations of the relative motions of the particular system 
under discussion, but they are not independent of each other. 
Indeed, if it were possible to overlook the conditional relations 
(6), instead of the (jp+i) constants which are sufficient for the 
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complete solution of the physical problem under consideration, 
we should have placed at our disposal Sp constants determinable 
from observation. 

These remarks are equally true when the method of the 
variation of parameters is employed. If a', n’ denote two variable 
parameters subject to the condition 


( 7 ) ... 

we shall have 

( 8 ) . . . 

or, 


/= 


7 V 3 

777 ' 


2tf/V . 


dn' 

~cLt 


+ 3»'V 2 


dcd_ 

dt 


= o, 


did _ 3 n r da' 
dt 2 a' dt 


as a necessary , but not a sufficient conditional relationship be¬ 
tween the parameters 7 , n' for the solution of the physical 
problem under discussion. It will also be necessary for the 
solution of the problem of relative motions of the (p+ i) material 
points that the constants introduced in the integration of the 
differential equations from which a' and n' are found should be 
such that they give the same numerical value for f when substi¬ 
tuted in the equation (7) for all the masses, v,, v 2 , &c. ; and it 
is only by showing that this condition is satisfied that the truth 
of Newton's law of universal gravitation can, strictly speaking, 
be established. 

The angular velocities 

?? l . n 2 . . . &c. 

and the linear quantities 

a x . a 2 . . . &c. 

can, therefore, only be accurately found from observation sub¬ 
ject to the conditions imposed by the identities (6), and the 
introduction of redundant constants into the theoretical expres¬ 
sions must be carefully avoided. 

But if the expressions for the coordinates 

x + 8x, y + z + 5 z, 

which are compared with the results of observation, are sensibly 
incomplete , or affected with sensibly erroneous empirical correc¬ 
tions ; and values of the angular velocities n l9 n 2i &c., and of the 
linear constants a u a 2 , (fee., or of their equivalents 
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are thus obtained from a discussion of a limited series of residual 
discordances by the method of least squares or otherwise, there 
is no mathematical necessity whatever, if Newton's law of gravi¬ 
tation is true, that the equations of condition (6) should be 
satisfied by the values of the constants thus found ; or any proof 
that the physical problem has been solved by the process followed, 
unless it can be shown that the values of f found from the dis¬ 
cussion of the residual discordances, in connection with the con¬ 
stants n ly n 2 , &c., V,, v 2 , &c., are sensibly the same for all the 
material points. 

But in the problems of the relative heliocentric motions of 
the centres of gravity of the planets, the motions under mathe¬ 
matical treatment are not directly dependent upon the rotation 
of the Earth about its axis, or on the linear dimensions of the 
Earth ; and it is possible in the discussions of these relative 
motions to clear the fundamental constant f from fallible errors 
of determination, and to render its value absolutely definite 
throughout our mathematical investigations by the selection of 
suitable units of length and of time. But when this course is 
followed, the angular velocity of the Earth about its axis and 
the linear dimensions of the Earth referred to the units selected 
to render f definite will have to be found from a discussion of 
observations, and the values which will be found from these dis¬ 
cussions will necessarily be liable to fallible errors of determina¬ 
tion. 

The difficulties thus presented are inherent in the physical 
problem under consideration, and they cannot be avoided by any 
mathematical artifices of change of constants or of variables. 
The close agreement between the values recently obtained for the 
constant of solar parallax from discussions of such very different 
geocentric motions as those of the Moon, Venus, Mars, and the 
asteroids, values which may be said to be confined within such 
limits as 8^’So and S'^o, is sufficient to prove that the outstand¬ 
ing errors of our geocentric coordinates are confined within small 
limits. But the exact values of the solar parallax cannot be 
obtained with certainty from such discussions, however accurately 
the geocentric distances of a planet from neighbouring stars or 
from the Sun are measured at any particular time, until the ratios 
of the linear constants 

C?2> &C., 

as well as the values of such constants as 

^ 2 , &C., 

are accurately known. I should be inclined to believe that 
the progress of gravitational astronomy is at present retarded 
from incompleteness of our theoretical investigations, the employ¬ 
ment of inconsistent conditions between the variable t and 
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facts of observation, and the introduction of redundant con¬ 
stants to diminish the residual errors formed from observations 
extending over only limited intervals of time, rather than from 
existing errors of observation. But such errors are undoubtedly 
quite sensible, and every step which tends to diminish them is 
one in the right direction, and should be most carefully utilised. 


Note on the Distribution of the Stars generally used for the 
Determination of Clock-Error. By E. J. Stone, M. A., F.R.S., 
Radcliffe Observer. 

The authorities at Greenwich for a great many years have 
annually published and distributed a list of the Mean Bight 
Ascensions of some 215 stars, brought up from the latest Green¬ 
wich Catalogue , and corrected, when necessary, from recent 
observations. This list has been found exceedingly useful by 
many astronomers, more particularly for the determination of 
the errors of their clocks. The apparent places of most of the 
Greenwich stars are included in the current Nautical Almanacs ; 
and I understand that the whole of them will be included in 
future volumes. But the distribution of the clock-stars in this 
list about the equator is not very satisfactory, and I have 
thought that this might easily be improved ; and, to test 
opinion on this point, the following list has been prepared. I 
have not ventured to cut out any of the stars which are now 
on the Greenwich list; but if the number of stars required to fill 
up the gaps should be thought to render the total number of 
stars inconveniently large for that constant re-observation which 
is desirable, there are some of the stars now in the list which 
might be cut out without any serious practical loss. The Bight 
Ascensions given for the additional stars are close approximations 
to the true values ; and although they are not at present as 
accurate as the Bight Ascensions of most of the other stars, the 
errors are so small that, if adopted for clock-stars, their places 
would be rapidly improved without leading to any practical in¬ 
convenience. 

Two diagrams showing the general distribution of the stars 
are appended. 
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